Exact Solution of an Octagonal Random Tiling Model 
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We consider the two-dimensional random tiling model introduced by Cockayne, i.e. the ensemble of all 
possible coverings of the plane without gaps or overlaps with squares and various hexagons. At the appropriate 
relative densities the correlations have eight-fold rotational symmetry. We reformulate the model in terms 
of a random tiling ensemble with identical rectangles and isosceles triangles. The partition function of this 
model can be calculated by diagonalizing a transfer matrix using the Bethe Ansatz (BA). The BA equations 
can be solved providing exact values of the entropy and elastic constants. 

PACS numbers: 61.44.-(-p, 05.20.-y, 64.70.Rh 



Since the discovery of quasicrystals, materials with 
non-crystallographic rotational symmetry and quasiperi- 
odic translational order have been modelled by tilings. 
A tiling model consists of a set of elementary building 
blocks, tiles, that cover space without gaps or overlaps. 
One of the main questions concerning quasicrystaUine al- 
loys is their thermodynamic stability. It has been argued 
by several authors [Qj^] that this may result partly from 
entropy associated with local random rearrangements of 
the tiles. One is then naturally led to study ensembles of 
'random tihngs' |Q. 

It is now known for some time [^[^] that in two di- 
mensions the square-triangle random tiling (RT) model, 
which has a twelve-fold rotational symmetry, can be 
solved, i.e. its entropy and phason elastic constants can 
be calculated exactly. In this letter we give the results of 
such a calculation for an eight-fold symmetric RT-model. 

The model under consideration consists of squares and 
hexagons of arbitrary size, and was first introduced by 
Cockayne The hexagons are built out of rectangles 
with sides 1 : V2 and a pair of isosceles and rectangular 
triangles. The squares can be viewed as two triangles. 
The model is therefore equivalent to a triangle-rectangle 
random tiling with an extra Boltzmann-weight such that 
the two ways two triangles form a square are counted as 

There is a repulsive potential of kT log 2 for each pair of 
triangles adjacent by their long edge. Since the configu- 
ration of vertices does not depend on which diagonal is 
drawn in a square, the partition sum is precisely the sum 
of all vertex configurations, rather than of all tilings. 

It is interesting to note that the perfect quasicrystaUine 
square-hexagon tiling generated by an inflation rule § is 
in one-to-one correspondence with the binary octagonal 
tiling of squares and rhombi. Although the random tiling 
ensemble of the latter set of tiles has been studied, d), 
no exact solution in the quasicrystaUine phase has been 
found yet. 



Like in the square triangle tiling we can set up a trans- 
fer matrix. This is done by decomposing the tiling into 
layers. Different layers are bounded by the short horizon- 
tal edges, the horizontal diagonals of the squares and the 
almost horizontal diagonals of the ±j tilted rectangles. 
In addition, the layer edges cut the triangles and rectan- 
gles with a vertical long edge in half. In this way the tiles 
are deformed in such a way that the vertices of the tiling 
form a subset of those of the square lattice, see Fig. |l|(b). 
The horizontal diagonals of the squares are denoted by 
the dashed lines in Fig. |l](a). A matrix element Tij of 
the transfer matrix T is if layer j can not be followed 
by layer i. Otherwise it is given by the statistical weight 
of the layer i. 

In the following we will call the tilted rectangles R± 
and the rectangles with the short and long horizontal 
edge by Rg and Ri respectively. 




FIG. 1. (a) Patch of the tihng model, (b) Corresponding 
patch on the lattice. Bold solid lines are domain walls of 
horizontal short edges, referred to as s- walls in the text. Bold 
dashed lines represent the Z- walls of horizontal long edges. 

Because different tiles of the original tiling are mapped 
onto the same shapes on the square lattice, we have to 
decorate the new configurations. This is done with bold 
dashed and solid lines, see Fig. |l|(b). Thus it is clear 
that the horizontal short and long edges of the original 
tiling form domain walls, which we denote by type s and 
I respectively. Between two layer edges on the square 
lattice, the s-walls step one unit to the left and the l- 
walls do not move. When two walls cross, the s-wall 



1 



may either jump over the Z-wall moving two places to 
the left and thereby creating a rectangle Ri, or, over two 
layers, the walls may exchange place creating a rectangle 
Rs ■ In the latter case the crossing therefore is completed 
after application of the transfermatrix twice. It may also 
happen that two walls of type s and one of type I cross 
simultaneously over two layers. The Z-wall and the s-wall 
nearest to it then exchange place, while the the second 
s-wall jumps over both these walls moving three places 
to the left creating a rectangle i?_ . 

We can express the tile densities in terms of the domain 
wall densities. We shall denote the horizontal size of 
the tiling by L and the corresponding system size of the 
lattice model by A''. Let A/s = i?; — i?^ = 0, i.e. both 
types of collisions of two domain walls occur with the 
same frequency, and let Us and rt; be the number of s- 
and Z-walls. Apply the transfer matrix p = 2N — Us times 
on some initial configuration of domain walls at t = 
on the lattice, and suppose that both types of collisions 
occur for every pair of s- and /-walls. The final state 
at t = 2N — ng will then be the same configuration of 
domain walls as the initial one shifted by n;. The total 
number of rectangles and of triangles per layer can then 
be calculated to be: 

nrect = N - Us - ni + 2nsni/p. 
ritri = 2{ns+ni) -6nsni/p. (1) 

The tile densities that belong to the quasicrystalline 
phase are Urect/N = 6 — 4-\/2, ntri/N = 12\/2 — 16, 
corresponding to an area fraction of triangles at = 1/2. 

As a function of the domain wall densities, the model 
displays two incommensurate phases. A 4-fold symmet- 
ric phase is formed in the high density region, at > 
1/2, where the triangles form octagonal and square cells 
bounded by domain walls consisting of rectangles. There 
is a 2-fold symmetric phase in the low density region 
where the rectangles from rectangular cells bounded by 
the domain walls consisting of triangles. 

Like the square-triangle tiling this tiling has the irro- 
tational property which implies that, according to the 
random tiling hypotheses, the entropy density (Ja = S/A 
has the following form: 

aa = (Tafl - (TrE)' + ]^K^ det E + O (E^) (2) 

where E is the phason strain tensor. The conditions on 
the elastic constants for E = to be a local maximum 
are: 

if^>0, K^>Q, AK^~K^>Q. (3) 

We denote the deviations of the ideal tile densities by 

5is = ni\p2- Ug. (4) 
A± = N -m-n,. (5) 



The quadratic forms in (g) can be expressed in these: 

(TrE)' = -1 (2<5;. - A± (2 - V2) - A,, (l + V2) ) ' . 

det E = (Sl - (2 - V2) A±^,, - 2A^V2 

- (1 + x/2) A,, [Su - (1 - I/V2) A±)) . 

(6) 

The quantities rtj and n; are conserved by the action of 
the transfer matrix T. To control the average value of 
A/s, the tiles Rs and Ri are given a weight exp(— ^) and 
exp((/)) respectively. Furthermore, as the tiles Rs and i?+ 
in the lattice representation have an area that is twice 
that of the other two transformed rectangles, we have to 
introduce a chemical potential for them to compensate 
for this asymmetry. The tiles -Rs and therefore get 
an extra weight exp(7y). 

The free energy per layer of the lattice model is given 
by the logarithm of the largest eigenvalue of T: 

F{ns,ni,(t)) = -log A 

= -S-<l>Ais 

-v{nR,+nR^). (7) 

We denote the horizontal coordinate of the ith s-wall 
by and of the kth /-wall by Zk- The vertical coordinate 
is denoted by by t. Let Z, be the total number of /-walls to 
the left of the ith s-wall, then the quantity £,i+t+li mod 2 
is conserved for every s-wall. This means that the s-walls 
lie on a sublattice structure and split up into two kinds: 
odd and even ones. Denote their coordinates by Xi and 
Tjj respectively. 

The eigenvectors of T function of the coordinates 
X, y and z of the domain walls are of the Bethe-Ansatz 
form. If the domain walls all are separated the Ansatz 
is: 

ir.^t.p i=l j=l k=l 

The form of the eigenvector for configurations where do- 
main walls cross can be found by application of T on (H) . 
Here, Wk = exp{iqk), Ui = exp(ipo^i) and Vj = exp(ipej) 
are the exponentiated momenta and p, tt and /i are the 
permutations of these belonging to the /- and odd and 
even s-walls respectively. The amplitudes A depend on 
the permutations p, tt and /i and on the configuration of 
the various domain walls. These together are coded in 
a vector F in the following way. Let r be the vector of 
coordinates Xi, yj and Zk of all domain walls, ordered so 
that r„i < rm+i- The entries of F are the elements of 
the permutations tt, and p. The order of succesion in 
F of elements taken from tt, and p matches that of the 
elements of x, y and z respectively in r. So, for example. 
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in the case of an odd s-wall at xi and an Z-wall at zi 
we would either have xi < zi, or xi > zi. In the first 
case we write r = (xi, zi) with T — (tti, pi), while in the 
second case r = {zi,xi) and F = (pi,7ri). 

When the different domain walls are separated, the 
transfer matrix shifts all s- walls to the left and leaves all 
Z-walls at rest, so the eigenvalue of T must be 



A 



n-n 



(9) 



Inspecting the eigenvalue equations for the case that an s- 
and an /-domain wall collide one sees that the amplitudes 
A before and after the collision must satisfy the following 
relation for (0) to be an eigenvector of T: 



Ai.. 



A{. ..pk^TTi . . .) 



(10) 



A same relation holds for the amplitudes with tt^ replaced 
by /ij and u replaced by v. From configurations involv- 
ing three domain walls one deduces that interchanging 
domain walls of the same kind in the amplitude gives a 
factor —1 and that interchanging an odd and an even 
s-wall leaves the amplitude unchanged. The eigenvalue 
equations therefore do not mix the momenta of the even 
and odd s-walls. It turns out that all relations among 
amplitudes involving more than two domain walls fac- 
torize into the ones already mentioned. These relations 
therefore suffice to make (||) an eigenvector of T. Im- 
posing periodic boundary conditions and eliminating the 
amplitudes A from the eigenvalue equations one gets the 
following equations for the momenta: 

uf = (-1)"=,-! Yl {e^u, + e^-^u^'w^') , 



k=l 
ni 



3 ""fc y ' 



fc=i 



W 



-Q (eV^ + e^-V^ii;^!) x 



1=1 



n 



(11) 



These are the so-called Bethe Ansatz equations (BAE). 
Like the BAE for the square-triangle tiling these equa- 
tions can be solved along a line in the thermodynamic 
limit for the largest eigenvalue. The details of this cal- 
culation, which resembles the one by Kalugin for the 
square-triangle tiling Q, will be published elsewhere. 
Here we only give the results. 

The entropy can be calculated exactly in the entire 
regime at > 1/2, A± = Ais = 0. In this region the 
three curves formed by the solutions of the BAE have 
the same limitpoint. After a change of variables, this 



limitpoint can be written in the notation of Kalugin as 
b = i\b\e^'^. The tile densities and the area fraction can 
be expressed in 7. 

n,/N = 1 - n^/N = ^-^'"^^'l (12) 



o-t 



1 + \/2cos7/2' 

/^^xV2_ccW2^ (13) 
V / 1 + cost/2 ^ ' 



The entropy per area of the square-hexagon random 
tiling in the regime 1/2 < at < 1 in terms of 7 is: 

a. ^ .,-^+^. r2^/21o, ' 



cos 7 



4 (cos 7/2 + 1) 

i /o • /oM /l + cos(7r/4 + 7/2)\ 

— cos 7/2 — sm7/2 log -. — J— 

^ " " ' s> \^i_cos(7r/4-F7/2) J 

1^ ■ /l + cos(7i-/4-7/2)\\ 

-(cos7/2-Hsin7/2 log ^4 

^ " " ' ^ Vl-cos(7r/4-7/2) / / 

(14) 

The entropy has its maximum at 7 = 0. Expanding CTq 
up to second order in 7 results in: 



O'a = CTa.fl — 1 



,l + \/2 
32V2 



(4-log4-y21og(l + %/2)) , 

(15) 



where the residual entropy at 7 = is given by 



2\/2 
0.1193642186 



(log4- V21og(l + \/2 



(16) 



As in the square triangle tiling, the entropy is a convex 
function of the area fraction at, see Fig. 0. 
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FIG. 2. (Ja as a function of ott. The solid line corresponds 
to the exact solution (p^. The dots are numerical results for 
N = 198. 
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In the 2-fold phase {at < 1/2) this exact calculation 
fails because the solution curves do not have the same 
limitpoint. Nevertheless, it is possible to calculate the 
lowest order correction to the entropy. With e ~ A±/N 
this is given by 

<Ja = <Ja,o - ^7^^ (4 - log4 - V2\og (1 + n/2)) 

-62i^flog4 + \/21og(l + V2)). (17) 



This expression gives the exact slope of the numerical 
curve shown in Fig. ^ for at t 1/2: 

^log4 + \/21ogfl + \/2)) . (18) 
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Expanding equation (|l^) up to first order in 7 and using 
(^ — (H), it is straightforward to find the elastic constants 
Kf, and from (H) and (|l|). Their numerical values 
are: 

Kf, = 0.2842712..., = 0.7366252.... (19) 

Since the two elastic constants fulfill the relations (|^) 
the quasiperiodic eightfold symmetric state is entropi- 
cally stable. 

In this paper we succesfuUy apply the Bethe Ansatz 
to an octagonal random tiling model. The BAE ( |ll| ) 
are solved to find exact values of the entropy and elastic 
constants. The model shows qualitatively the same be- 
haviour as the square-triangle tiling. It is not yet clear 
how generic the solvability of these two tilings is, but we 
have discovered that a 10-fold symmetric tiling of rectan- 
gles and triangles does admit a Bethe Ansatz. It appears 
however from numerical calculations that their solutions 
do not allow for an exact solution using the method of 
Kalugin which is employed in this paper. 
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